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Abstract 

Background: SKIP is a transcription cofactor in many eukaryotes. It can regulate plant stress tolerance in rice and 
Arabidopsis. But the homolog of SKIP protein in soybean has been not reported up to now. 

Results: In this study, the expression patterns of soybean GAMYB binding protein gene {GmGBPl) encoding a 
homolog of SKIP protein were analyzed in soybean under abiotic stresses and different day lengths. The expression 
of GmGBPl was induced by polyethyleneglycol 6000, NaCI, gibberellin, abscisic acid and heat stress. GmGBPl had 
transcriptional activity in C-terminal. GmGBPl could interact with R2R3 domain of GmGAMYBl in SKIP domain to 
take part in gibberellin flowering pathway. In long-day (16 h-light) condition, transgenic Arabidopsis with the 
ectopic overexpression of GmGBPl exhibited earlier flowering and less number of rosette leaves; Suppression of 
AtSKIP in Arabidopsis resulted in growth arrest, flowering delay and down-regulation of many flowering-related 
genes {CONSTANS, FLOWERING LOCUS T, LEAFY); Arabidopsis myb33 mutant plants with ectopic overexpression of 
GmGBPl showed the same flowering phenotype with wild type. In short-day (8 h-light) condition, transgenic 
Arabidopsis plants with GmGBPl flowered later and showed a higher level of FLOWERING LOCUS C compared with 
wild type. When treated with abiotic stresses, transgenic Arabidopsis with the ectopic overexpression of GmGBPl 
enhanced the tolerances to heat and drought stresses but reduced the tolerance to high salinity, and affected the 
expressions of several stress-related genes. 

Conclusions: In Arabidopsis, GmGBPl might positively regulate the flowering time by affecting CONSTANS, 
FLOWERING LOCUS T LEAFY and GAMYB directly or indirectly in photoperiodic and gibberellin pathways in LDs, but 
GmGBPl might represse flowering by affecting FLOWERING LOCUS C and SHORT VEGETATIVE PHASE in autonomous 
pathway in SDs. GmGBPl might regulate the activity of ROS-eliminating to improve the resistance to heat and 
drought but reduce the high-salinity tolerance. 
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Background 

Flowering time plays key roles in plant development, plant 
adaptation to growing regions [1], crop yield [2,3] and dis- 
ease resistance [4]. Abiotic stress such as drought, high sal- 
inity and heat, is another major limiting factor for crop 
yield [2]. Genetic and molecular analyses has revealed that 
several distinct but linked signaling pathways regulate the 
flowering-time in response to light and temperature signals 
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or to internal signals such as vernalization, autonomous 
pathways and gibberellin signal pathway [1,5-7]. Interaction 
of different signal transduction pathways before their con- 
vergence may allow a coordinated regulation of the activity 
of the respective pathways. The signals from these pathways 
are integrated through up-regulating the expression of one 
or more common target genes: FT, SOC1 and LFY to deter- 
mine when flowering. The expressions of SOC1 and FT are 
negatively regulated by FLC but positively regulated by CO 
[8,9]. FLC, a MADS -box transcription factor, is a negative 
regulator of floral initiation and an integrator of the 
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autonomous and vernalization pathways. CO, encoding a 
B-box zinc-finger transcription factor, integrate circadian 
clock and light signals and up-regulate the expression of 
FT and SOC1 directly in the long-day pathway [10]. 
CONSTANS and GAMYB (AtMYB33) regulate the expres- 
sion of LFY by binding with different cis-acting elements in 
the promoter of LFY respectively [11-14]. SKIP is a tran- 
scription cofactor in many eukaryotes. All the SKIP homo- 
logs identified so far contained a SKIP domain with an 
S-N-W-K-N peptide signature and may have conserved 
basic functions, such as acting as a cofactor in transcrip- 
tion and splicing [15]. However, the derived or additional 
functions of the SKIP homologs varied among species. 
SKIP was an essential protein for pre-mRNA splicing in 
Saccharomyces cerevisiae [16]. In Drosophila melanogasten 
SKIP was involved in ecdysone-stimulated transcription 
[17] and known to be a coactivator in Notch [18,19]. In 
Caenorhabditis elegans, SKIP was an essential component 
of many RNA polymerase II transcription complexes and 
indispensable for C. elegans development [20]. In Hordeum 
vulgare L., SKIP could interact with GAMYB [21]. In Oryza 
sativa L., OsSKIPa could positively modulate cell viability 
and stress tolerance [22]. AtSKIP functions as not only a 
positive regulator and putative potential transcription factor 
in the abiotic stress signaling pathway in Arabidopsis [23], 
but also a component of the spliceosome linking alternative 
splicing and the circadian clock in Arabidopsis [24]. How- 
ever, there has been no report on the identification of SKIP 
homologs in soybean and on their functions in the regula- 
tion of flowering time in any species so far. 

In this study, the expression patterns of soybean 
GAMYB binding protein encoding a homolog of SKIP 
protein were analyzed under several abiotic stresses and 
light conditions. Yeast two-hybrid assay was performed 
to identify the interactions between GmGBPl and 
GmGAMYBl. The phenotypes of transgenic Arabidopsis 
lines in different stress and day-length were analyzed to 
study the function of GmGBPl. 

Results 

GmGBPl encodes a homolog of the human SKIP 
transcriptional coregulator 

The SKIP homolog gene in soybean was obtained from 
subtracting long-day from short-day treated mRNA. 
The predicted gene, named as GmGBPl (Glycine max 
GAMYB-b'm&mg protein gene, GenBank DQ1 12540), 
wasamplified from the cDNA of DongNong 42 using 
PCR. The full length cDNA of GmGBPl, containing 
2253 bp with an open reading frame of 1,839 bp, was 
predicted to encode 612 amino acids. GmGBPl protein 
had a predicted isoelectric point (pi) of 8.69 and a mo- 
lecular weight (MW) of 69.10 kDa. 

Sequence alignment indicated that SKIP proteins were 
highly conserved in eukaryotes. The SKIP protein domain 



of GmGBPl was located between amino acids 190 and 
360 (Figure 1A) with an S-N-W-K-N peptide signature 
(Figure IB). A BLAST search was performed against 
the assembled soybean genome databases (http://www. 
phytozome.org) to analyze the SKIP sequences in viridi- 
plantae. The phylogenetic tree based on the SKIP sequences 
agreed well with the evolutionary relation among these spe- 
cies those were divided into four groups (Figure ID). 

A yeast transformation was used to detect the tran- 
scriptional activity of GmGBPl. GmGBPl, GmGBPla, 
GmGBPlb and GmGBPl c were fused to the GAL4 
DNA-binding domain and transformed into yeast strain 
YRG-2 cells respectively. The proteins of GmGBPl and 
GmGBPl c were both capable of inducing LacZ expres- 
sion in yeast cells and could grow in the medium lacking 
His, indicating that GmGBPl protein had transcriptional 
activity in its C-terminal domain (Figure 1C). 

GmGBPl expression in soybean 

A quantitative RT-PCR was performed to analyze the ex- 
pression of GmGBPl in soybean. GmGBPl could express 
in all tissues or organs investigated in both LDs and SDs, 
although the expression levels were different. The expres- 
sion levels in stem and immature seed were much higher 
than those in root, trifoliate leaf, flower bud, pod and leaf. 
The trifoliate leaf had the lowest expression level in LDs, 
but the highest expression level in SDs in the seven sam- 
ples detected. All the organs, except of stem, showed 
higher expression in SDs than in LDs (Figure 2F). The ex- 
pression of SKIP homologs in rice and Arabidopsis have 
been confirmed to be affected by stresses and phytohor- 
mones (Xin et al. 2009; Gah-Hyun et al. 2010). Therefore, 
the expression level of GmGBPl in soybean was also ana- 
lyzed with stresses and phytohormones treatments in the 
present study. The results indicated that the transcrip- 
tional level of GmGBPl was increased after treated with 
100 umol ABA (Figure 2A), 100 umol GA (Figure 2C), 
200 mM NaCl (Figure 2B), 8% PEG6000 (Figure 2D) and 
temperature 42°C respectively (Figure 2E). 

GmGBPl interacts with GmGAMYBl in yeast cells 

In barley, it was confirmed that HvGAMYB could inter- 
act with a SKIP domain by yeast two-hybrid system 
[21]. Therefore, a yeast two-hybrid assay was performed 
to study the interaction between GmGBPl and 
GmGAMYBl. GmGAMYBl contained a typical R2/R3- 
MYB DNA-binding domain and three conserved regions 
(BOX1, BOX2 and BOX3). GmGAMYBl was more 
homologous to AtMYB33 rather than AtMYB65 and 
AtMYBlOl (Figure 3A, Figure 3B). As GmGBPl showed 
transcriptional activity in the C-terminal domain, the 
pBD-GmGBPlz and pBD- GmGBPlb were chosen for 
the bait vectors to search the protein interaction 
with GmGAMYBl. Moreover, the prey was divided 
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Figure 1 Identification of SKIP homologs in soybean. (A) Predicted domains and motifs of soybean SKIP homolog, GmGBPl. (B) Multiple 
alignment of the predicted amino acid sequences in peptide signature of SKIP homologs from soybean {GmGBPl, DQ1 12540), Arabidopsis (AtSKIP, 
O80653), rice {OsSKIP), human (HsSKIP, Q13573), Coenorhoobditis elegons (CeSKIP, Q22836), Drosophilo (DmBX42, P39736), Socchoromyces cerevisioe 
{ScPRP45, Q09882). (C) Detection for self-transactivation of GmGBPl. The deletion mutants (shown on right), GmGBPl, GmGBPl a (amino acids 1- 
189), GmGBPIb (amino acids 190-356) and GmGBPl c (amino acids 357-612), were constructed by PCR. Self-transactivation, as illustrated by the 
color reaction in the 5-bromo-4-chloro-3-indolyl-(3-D-galactopyranoside(X-Gal) assay, occurred with In for pBD-GmGBPl and pBD-GmG£P7c 
constructs. Overnight incubation failed to produce a self-activation reaction with the pBD-GmGBPla and pBD-GmGBPlb constructs. And the 
transformants were streaked on SD agar plates lacking His, Trp and incubated at 30°C for 72 h. The pBD-GmGBPl and pBD-GmGBPlc constructs 
could grow normally but no growth for the pBD-GmGBPIa and pBD-G/nG£P7b constructs. (D) Phylogenetic relationships of SKIP homologs in 
plant constructed using the neighbor-joining method with the program MEGA 5.05. 



into two parts (1-141,142-538) by the domain area 
of R2R3, as pAD-GmGAMYBl, pAD-GmGAMYBla and 
pKD-GmGAMYBlb. The co-transformation with pBD- 
GrnGBPla and the prey plasmid could not grow on SD/- 
TLH, indicating that the section of GmGBPl had no ability 
to interact with GmGAMYBl (Figure 3C). Furthermore, 
the yeast cells with pBD-GmGBPlb and pAD-GmGA- 
MYBlb showed no ability to interact with GmGAMYBl 
(Figure 3D). In contrast, the section b of GmGBPl inter- 



acted with GmGAMYBl by binding to the domain of 
GmGAMYBl as the result of the growth on SD/-TLH, 
exhibiting the positive color reaction (blue on the X-Gal 
filter assay, Figure 3B). 

The homozygous mutant atskip of Arabidopsis was lethal 
at late stage 

An insertion of a T14N5 fragment at the 1344 bp of 
exon of AtSKIP caused mutant (Figure 4A). The mutant 
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Figure 2 Expression patterns of GmGBPI in soybean. (A, B, C, D and E) Expression level of GmGBPl in leaf of soybean treated with 100umol 
ABA, 200mM NaCI, 100umol GA, 8% PEG6000 and heat stress respectively. (F) Tissue-specific expression of GmGBPI in LDs and SDs. Tissues tested 
were root (RO), stem (ST), leaf (LE), trifoliate leaf (TL), flower bud (FB), pod (PO) and immature seed (IS). 
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atskip had the same germination rate with the WT plant 
(>99%), and all of the plants could grow normally in the 
MS medium. However, 100% of the homologous plants 
died within 2 weeks after transplanting in soil before flow- 
ering whereas all WT and heterozygous plants grew nor- 
mal under the same conditions (Figure 4B, Figure 4C). The 
homozygosity of mutant was detected by PCR and RT- 
PCR after transplanting in soil. 

Considering the homozygous mutant of AtSKIP was le- 
thal, RNAi was used to produce the atskip-i plant to sup- 
press the expression of this gene. Although the expression 
of AtSKIP was only partially suppressed (Figure 4E), the 
plants showed abnormal phenotypes such as late flower- 
ing, small rosette leaves. According to RT-PCR analyses, 
the expressions of several key flowering- related genes were 
down-regulated (Figure 4E). The promotion and inhibitory 
factors of flowering were both down-regulated, which sug- 
gested that SKIP might be involved in more than two dif- 
ferent flowering pathways. 

Ectopic expression of GmGBPI alters stress tolerance in 
Arabidopsis 

Since GmGBPI was up-regulated by stresses and stress- 
related phytohormones, Arabidopsis thaliana plants, con- 
taining extra copies of GmGBPI gene, was produced to 



investigate whether overexpression of this gene could en- 
hance stress tolerance. The transgenic lines were con- 
firmed by PCR and RT-PCR with homolog-specific 
primers (data not shown). In the MS medium containing 
ABA (1 and 3 |iM), the germination rates of GmGBPI -ox 
and atskip-i plants were both reduced compared to 
WT (Figure 5A), indicating that ectopic expression of 
GmGBPI could increase the growth-suppression effect of 
ABA and AtSKIP might have the ability of reducing the 
negative effect of ABA. In the MS medium containing 
NaCI (150 and 200 mM), the germination rates of 
GmGBPl-ox and atskip A plants were also lower than the 
WT plants (Figure 5A, Figure 5C), suggesting that 
GmGBPI decreased salt tolerance and AtSKIP improved 
the ability to salt tolerance. However, the germination 
rates of the three plants under 300 mM mannitol showed 
no difference (data not shown). A further study was per- 
formed to detect the stress tolerance of GmGBPI. When 
20 days old plants growing in soil under SDs were treated 
with drought condition (no watering for 7 days), nearly all 
of the GmGBPI -ox plants were recovered at 3 days after 
re-watering, whereas only 40-60% of WT plants and none 
of atskip A plants were recovered (Figure 5B). The root 
lengths of both GmGBPl-ox and atskip-i plants were 
shorter than that of WT plants in the MS medium treated 
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Figure 3 Yeast Two-Hybrid Assay of Interaction between GmGBPI and GmGAMYBI. (A) Phylogenetic relationships between GmGAMYBI 
(HM447241) and GAMYB homologs (AtMYB33, NMJ80448; AtMYB65, NM_1 1 1977; AtMYBWl, NM_1 28805.3) in Arabidopsis using the neighbor- 
joining method with the program MEGA 5.05. (B) Multiple alignment of the amino acid sequences in domain and motify of GmGAMYBI , 
AtMYB33, AXMYB65 and AtMYBWl. (C) The growth of Yeast YRG-2 cells cotransformated by pBD-GmGBP1a with pAD-GmGAMYBl , pAD- 
GmGAMYBIa (amino acids 1-141) or pAD-GmGAMYBlb (amino acids 142-538) on synthetic minimal (SD) medium lack of Trp and Leu (SD/-TL) or 
Trp, Leu and His (SD/-TLH). (D) The growth of Yeast YRG-2 cells cotransformated by pBD-GmGBPlb with pAD-GmGAMYBl , pAD-GmGAMYB1a or 
pAD-GmGAMYBlb on SD/-TL and SD/-TLH, and the X-Gal assay for the transformations. 



with 150 mM NaCl for 10 days (Figure 5D); The 
GmGBPI -ox and atskip-i plants treated with 200 mM 
NaCl for 10 days under LDs were nearly dead, while 
80-90% of WT plants survived (Figure 5D). 

Heat tolerance experiment revealed that GmGBPI 
might improve heat tolerance of Arabidopsis plants 
with ectopic expression of GmGBPI. 7days old plants 
of transgenic Arabidopsis, treated with 37°C for 5 days, 
exhibited symptom of flagged leaves, but WT and 
atskip-i plants showed more serious symptom with 
whitened leaves (Figure 5E: a, b and c). Furthermore, 
WT and atskip-i plants of 15days old showed more 
whitened leaves, even dead after treated with 37°C for 
8 days, whereas GmGBPI -ox plants showed more vital 
(Figure 5E: d, e and f). 



The GmGBPI -ox plants had a higher Superoxide dis- 
mutase (SOD) level than both atskip-i and WT plants 
under heat or drought conditions. The malondialdehyde 
(MDA) level of GmGBPI -ox plants was much lower 
than both atskip-i and WT plants, indicating that the ec- 
topic expression of GmGBPI in Arabidopsis might en- 
hance the heat and drought tolerance. However, the 
highest MDA level and the lowest SOD level of 
GmGBPI -ox plants implied the reduced salt tolerance in 
GmGBPI -transgenic Arabidopsis (Figure 5F). 

Stress resistance-related genes are affected by GmGBPI 

To gain further insight into the mechanism of the 
altered stress resistance of GmGBPI -ox plants, transcript 
levels of 8 stress-responsive genes [25-28] were assayed 
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Figure 4 Identification of atskip mutant plants in phenotypes and genotype. (A) A T-DNA insertion caused the mutant of AtS KIP. (B) 
Phenotypes of homozygous atskip mutant (HM), heterozygous atskip mutant (HZ) and wild type Col-0 plant (WT). Bars= 1.0cm. (C) PCR tests for 
genotype to identify the mutant. (D) RT-PCR analysis for the down-regulated genes {AtFT, NM_1 05222.2; AtCO, NM_001 03681 0.1; AtLFY, 
NM_1 25579; AtFKF, Q9C9W9; AtSOC 1, NP_1 82090.1 ; AtFLC, Q9S7Q7; AtSVP, Q9FUC1) in atskip-i plants. 



in the GmGBPl -ox, atskip A and WT plants under normal 
and stress conditions. The expression levels of all the 8 
genes in GrnGBPl-ox were higher than WT and atskip-i 
plants (Figure 6). Moreover, ZAT7, ZAT12, HSFAlb and 
HSFA2 showed significantly higher expression level in 
GrnGBPl-ox than in WT and atskip-i plants after treated 
with 42C for 2h (Figure 6), MYB2, MYB96, MYC2, RD26 
and ZAT12 were up-regulated after drought for 7days 
(Figure 6). These results suggested that overexpression of 
GmGBPl can increase the transcript levels of some stress- 
resistance genes under drought or heat stress conditions, 
thus improving the drought or heat resistance of the 
transgenic plants. However, the expression levels of RD26, 
ZAT7 and ZAT12 in GmGBPl -ox plants were significantly 
lower than WT and atskip-i plants after treated with 200 
mM NaCl for 7 days (Figure 6), indicating that GmGBPl 
might reduce the resistance to high salinity. 

Early flowering of GmGBPI-ox plants in long-day 
condition 

The response to day-length implied the possibility of 
flowering-related functions of GmGBPl. In LDs condition 
(16-h light), the T3 plants of GmGBPI-ox Arabidopsis 
flowered 2 days earlier and had fewer rosette leaves than 
WT plants, and the atskip-i plants flowered significantly 
late (Figure 7 A, Figure 7D). Furthermore, the expressions 
of CONSTANS (CO) and FLOWERING LOCUS T (FT) in 
GmGBPI-ox plants were increased compared to WT and 
atskip-i plants, and atskip-i plants had the lowest expres- 
sion levels of CO and FT (Figure 7B), indicating that 



GmGBPl could control the flowering time by up-regulating 
flowering-related gene CO and FT in LDs. 

However, in the short-day (SD) condition, the T3 plants 
of GmGBPI-ox Arabidopsis flowered latest and had the 
most rosette leaves (Figure 7A, Figure 7D). In contrast, 
atskip-i plants showed the earliest flowering and had the 
fewest rosette leaves. This result suggested that there 
might be some other ways to regulate flowering in SDs. 
The key flowering inhibitor gene, FLC, was significantly 
up-regulated in GmGBPI-ox plants (Figure 7B). All the 
results suggested autonomous pathway could account for 
the late flowering phenotype. 

GAMYB genetically acts downstream of GBP to regulate 
flowering time 

Previous reports showed that AtMYB33 influenced flower- 
ing by mediating GA responsiveness of the LFY promoter 
[12,14], and was down-regulated in atskip-i plant in this 
study (Figure 4D). myb33 mutant plants (SALK058312) 
and myb33 mutant plants containing extra copies of 
GmGBPl gene were introduced to investigate the impact of 
the interaction between GmGBPl and GmGAMYBl on 
flowering. In LDs, no effect on flowering time and leaf 
numbers was observed in GmGBPI-ox/ my b33 plants 
(Figure 7A, Figure 7D), and the expression level of LFY in 
GmGBPI-ox/ myb33 plants was same to the level of 
WT, but GmGBPI-ox plants had a high level of LFY 
(Figure 7C). GmGBPl was directly located on the upstream 
of GAMYB, and regulate flowering in LDs. 
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Figure 5 Performance of seedling under stress condition. (A) The germination rate of plants treated with 1 \M ABA, 3 uM ABA, 150 mM 
NaCI or 200 mM NaCI individually. (B) Performance of seedlings treated with drought. Effect of 20 days old plants growing in soil in SDs treated 
with drought (no water for 7days), photos at 3 days after rewatering, Bars= 1.0 cm. (C) Phenotypes of seedlings treated with 150mM NaCI. 
Germination of plants growing in MS medium with 150mM NaCI, photos at 7 days after germination; Growth of plants cultivated in MS medium 
with 150mM NaCI vertically, photos at 10 days after germination. Bars= 1.0 cm. (D) Effect of 20 days old plants growing in soil in LDs treated with 
200mM NaCI every other day for 3 times, photos at 10 days after treated. Bars = 1.0 cm. (E) Performance of seedling under heat: a, b and c: effect 
of plants of 7 days old under 37°C for 5 days; d, e and f: effect of plants of 15 days old under 37°C for 8 days. Bars = 0.5 cm. (F) The relative level 
of MDA and SOD of plants under stress conditions. 20 days old plants growing in soil in SDs treated with 42°C for 36h, samples were prepared 
from rosette leaf after recovered for 1 day; 20 days old plants growing in soil in LDs, samples were prepared from rosette leaf after treated with 
drought (no water for 7 days); 20 days old plants growing in soil in LDs treated with 200mM NaCI every other day for 3 times, samples were 
prepared from rosette leaf after 7days. Data were from more than 8 individuals for each genotype ± SD. * f test, P < 0.05; ** t test, P < 0.01 . 



Since GmGBPl regulated flowering in autonomous 
pathway in SDs, a further study was performed to inves- 
tigate the responsibility of MYB33 in the late flowering 
phenotype. The result revealed that GmGBPl -ox/ myb33 
plants also showed late flowering like GmGBPl-ox 
plants (Figure 7 A and Figure 7D), suggesting that 
GmGBPl was independent of MYB33 in delaying flower- 
ing time in SDs. 



Discussion 

Essential roles of SKIP for plant growth and stress 
tolerance 

SKIP protein, as a spliceosome component, played a vital 
role in maintaining cell viability in yeast (PRP45), C. elegans 
(CeSKIP) and O. sativa (OsSKIPa) [16,20,22]. The homozy- 
gous mutant of AtSKIP in Arabidopsis resulted in death of 
plants, which was similar to the phenotype of severe growth 
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Figure 6 Expression of stress resistance-related genes in WT, GmGBPI-ox and atskip-'i plants. ZAT7, NM_1 14478.3; ZAT12, AY050915.1; 
HSFAlb, NMJ21688; HSFA2, NM_001 124916; RD26, NM_001 084983; MYC2, NM_102998; MYB96, NMJ25641; MYB2, NM_1 30287. 



arrest and even death of O. sativa with suppression of 
OsSKIPa and embryonic arrest of the CeSKIP mutant in 
C. elegans, which indicated that AtSKIP was also required 
for maintaining cell viability and normal growth in Arabi- 
dopsis. The indispensable role of SKIP homolog in keeping 
normal cell viability and growth might be conserved in 
plants, which was supported by that the plant SKIP homo- 
logs possessed the conserved SKIP domain for cell viability 
identified in PRP45. However, it needs further study to 
confirm whether the SKIP homologs in soybean could 
have the conserved function for cell viability. 

Drought and high-salinity repressed plant growth and 
limited seed yield. SKIP homologs in Arabidopsis and rice 
could improve the tolerance to drought and high-salinity 
[22,23]. In this study, the ectopic expression of GmGBPl in 
Arabidopsis could enhance the tolerance to drought. How- 
ever, GmGBPl-ox plants showed a salt-sensitive appear- 
ance, although atskip-i plants reduced the tolerance to salt 
as previous report [23]. Thus, the tolerance of SKIP homo- 
logs to salt might be not conserved in plant. Heat, one of 



serious environmental stresses, affected the growth of plants 
and the productivity of crops. There was no report about 
heat tolerance of SKIP homologs previously. Our result 
showed that GrnGBPl-ox plants increased the tolerance to 
heat whereas atskip-i plants reduced the tolerance to heat, 
indicating that AtSKIP has the ability of heat tolerance. 

Reactive oxygen species (ROS) could be induced by abi- 
otic stresses, and over-accumulation of ROS could lead to 
cell damage and even death. SOD was used to eliminate 
ROS, and MDA was the intermediate product during the 
elimination of ROS [29,30]. In this study, the raised con- 
tent of ROS in GmGBPl -ox plants treated with salinity 
and the decreased content of ROS in GrnGBPl-ox plants 
treated with heat (or drought) indicated that the altered 
stress tolerance of GmGBPl -ox plants may be partially 
due to the regulation of the activity of ROS -eliminating. 

GmGBPl regulates flowering time 

Interestingly, in this study ectopic expression of GmGBPl 
in Arabidopsis induced earlier flowering in LDs and late 
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D 

LDs SDs 



Genotype 


Days 


No. of 


Days 


No. of 




to Flower 


Rosette Leaf 


to Flower 


Rosette Leaf 


Wild type 


22.15 ±0.69 


9 ± 0.32 


38.96 ± 3.51 


16.42 ± 1.24 


myb33 


22.18 ±0.32 


9.3 ±0.24 


39.23 ± 2.92 


16.82 ± 0.96 


GmGBPIoxl 


20.23 ± 0.86* 


7 ± 0.44* 


49.34 ±6.13" 


26.64 ± 2.46" 


GmGBP1-ox2 


19.26 ±0.59' 


6.2 ± 0.64" 


46.22 ± 4.34" 


24.97 ± 1.96" 


atskip-rt 


26.17 ± 2.14* 


12 ± 0.78" 


33.82 ± 3.23" 


12.28 ± 1.32* 


atskip-\2 


28.18 ± 1.68" 


12.5 ± 0.86* 


32.12 ±2.24" 


11.93 ± 0.95" 


Gm GBP1-oxlmyb33A 


22.36 ± 0.96 


9.12 ±0.46 


48.26± 2.98" 


26.92 ± 1.54" 


Gm GBP1-ox/myb33-2 


22.29 ± 0.46 


8.94 ±0.72 


47.42 ±3.15" 


25.12 ± 0.93" 



Figure 7 Effect of WT, atskip-i, GmGBPI-ox, myb33 and GmGBP1-ox/myb33 plants. (A) Phenotypes of plants in LDs and SDs. Photos at 25 
days after germination in LDs; Photos at 42 days after germination in SDs. Bars= 1cm. (B) Expression level of flowering related gene in 
Arabidopsis. RT-PCR analysis for AtCO and AtFT in WT, atskip-\, GmGBPI-ox plants in LDs; RT-PCR analysis for AtFLC for plants in SDs. Samples were 
prepared from rosette leaves before bolting in LDs or SDs. (C) Expression level of AtLFY in plants. Samples were prepared from rosette leaves 
before bolting in LDs. (D) Flowering time of plants in LDs and SDs. Flowering time was measured as the total number of rosette leaves or days to 
bolting. Data were from more than 20 individuals for each genotype ±SD. WT plants compared with atskip-i, GmGBPI-ox, myb33 and GmGBPl- 
ox/myb33 plants, respectively. * showed significant to WT plant at 0.05 level (P < 0.05) by t-test; ** showed significant to WT plant at 0.01 level (P 
< 0.01) by t-test. 



flowering in SDs. The finding of a flowering-regulation 
function from a SKIP homolog has not yet been reported 
in any plants so far. The specific function of GmGBPl in 
regulating flowering time might be especially useful for 
developing soybean cultivars with the adaptability to broad 
grow regions. 



Flowering time was controlled by several signaling 
pathways, such as the day-length, vernalization, autono- 
mous pathways and gibberellin signal pathway [1,5-7], 
and it could be measured by scoring the number of ros- 
ette leaves at flowering time and the number of days 
from germination to bolting in Arabidopsis [31]. The 
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number of rosette leaves at flowering time and the num- 
ber of days from germination to bolting measured in 
GmGBPl -ox plants were both reduced in LDs, suggesting 
that ectopic expression of GmGBPl induced early flower- 
ing, atskip-i plants delayed flowering in LDs supported the 
early flowering of GmGBPl -ox plants in another hand. 
However, when GmGBPl ectopicly expressed in myb33 
plants, no early flowering could be observed in LDs, indi- 
cating that the early flowering function of GmGBPl might 
depend on the existence of MYB33. When the transgenic 
lines were transferred into SDs, the phenotypes of flower- 
ing were all changed. The flowering times of GmGBPl -ox 
and GmGBPl-oxlmyb33 plants were both delayed signifi- 
cantly, and atskip-i displayed the early flowering pheno- 
type in SDs. 

The expression levels of numerous flowering-related 
genes could be induced by day-length or gibberellin. The 
expression level of GmGBPl was regulated by both day- 
length and gibberellin, suggesting that GmGBPl might 
participate in both day-length and gibberellin signal 
pathways. In atskip-i plants of Arabidopsis, the expres- 
sion levels of more than 10 flowering-related genes were 
affected. In particular, the mRNA levels of several flow- 
ering integrators (FT,LFY and SOC1) were significantly 
lower. Our results showed that FT and LFY were both 
up-regulated in GmGBPl -ox plants in LDs, but LFY 
showed no change when GmGBPl was ectopic expres- 
sion in myb33 mutants. Considering the interaction be- 
tween GmGBPl and GmGAMYBl, and the conserved 
function on photoperiodic flowering of FT homologs in 
soybean [32-34], GmGBPl might regulate flowering time 
in two ways by its function as a transcription factor and 
interaction in LDs. On one hand, GmGBPl could regu- 
late the expressions of flowering-related genes in day- 
length signal pathway, and on the other hand, GmGBPl 
could bind with GmGAMYBl in gibberellin signal path- 
way to control flowering time. 

Many reports had revealed that GA pathway played a 
key role in flowering under SD condition when other 
regulatory pathways that promoted flowering were 
not active [1,7,35-37]. However, the late flowering of 
GmGBPl -ox plants indicated that there might be an- 
other way to control flowering in SDs. FLC was a nega- 
tive regulator of floral initiation and an integrator of the 
autonomous and vernalization pathways. FLC could dir- 
ectly down- regulated FT and SOC1 to repress the flower- 
ing in these pathways [8,9,38]. SVP, a MADS box 
transcription factor, could interact with FLC and acted 
as partially redundant repressors of flowering time with 
FLC [39]. The negative action on the phenotypes of the 
down-regulated FLC and SVP in atskip-i plants in LDs 
suggested that flowering time was regulated crossly by 
several signal pathways and SKIP might take part in 
more than two ways. In SDs, the expression of FLC was 



up-regulated in GmGBPl -ox plants, and was down- 
regulated in atskip-i plants as that in LDs. The late flow- 
ering of GmGBPl -ox plants in SDs elucidated that 
GmGBPl might delay flowering time in SDs through au- 
tonomous pathway by improving the expression level of 
FLC, a key flowering inhibitor factor. 

In general, GmGBPl might regulate flowering time by 
three signal pathways. GmGBPl positively controlled the 
flowering time by regulating CO, FT, LFY and GAMYB 
directly or indirectly in photoperiodic and gibberellin 
pathways in LDs, GmGBPl repressed flowering by regu- 
lating FLC and SVP in autonomous pathway in SDs. 

Diverse functions of SKIP homologs in plant 

SKIP had a conserved SKIP domain with an S-N-W-K-N 
peptide signature, and was considered as a cofactor for 
transcription regulation in all eukaryotes so far. How- 
ever, the derived or additional functions of the SKIP 
homologs varied among species. Transgenic rice that 
overexpressed OsSKIPa exhibited stress tolerances (abscisic 
acid, salt, mannitol) at both seedling and reproductive 
stages [22]. Overexpression of the AtSKIP gene in Arabi- 
dopsis modulated the induction of salt tolerance, dehydra- 
tion resistance and insensitivity towards abscisic acid under 
stress conditions [23]; However, ectopic expression of 
GmGBPl in Arabidopsis reduced tolerance to NaCl, but 
increased tolerance to drought and heat in this study. 

Suppression of OsSKIP resulted in growth arrest of rice 
due to the reduced cell viability in the active growth 
regions [22]. A decrease in AtSKIP expression led to 
altered plant development with the phenotype of reduced 
inflorescence stems and smaller rosette leaves [23]. In this 
study, the homozygous atskip mutant was lethal to the 
growth, and the knockdown of AtSKIP showed late flower- 
ing and increased number of rosette leaves. The ectopic 
expression of GmGBPl induced premature flowering of 
Arabidopsis in LDs. All the study on SKIP homologs in 
plant revealed that SKIPs might play a vital role in the 
growth and development of plant, but have different func- 
tion on the stress tolerance and flowering in plant. 

Plant SKIPs were divided into four groups based on 
the SKIP sequences by phylogenetic tree (Figure ID). 
The different groups among OsSKIP, AtSKIP and 
GmGBPl indicated the diversity of their functions. SKIP 
could be induced by various abiotic stresses, phytohor- 
mones treatments and day-length, but the expression 
patterns of SKIP varied among rice, Arabidopsis and 
soybean [22,23]. The transcriptional levels of OsSKIP, 
AtSKIP and GmGBPl were all up-regulated by ABA, 
NaCl and drought (PEG6000 or mannitol), but had dif- 
ferent scales. GmGBPl could be induced by day-length 
and heat, but no related report for OsSKIP and AtSKIP 
in rice and Arabidopsis. The various functions of SKIP 
homologs might be also owing to the diversification of 
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57<7P-interacting proteins SIPs. Previous studies showed 
that the SKIP homologs, PRP45 (yeast), BX42 (Drosophila), 
CeSKIP (C elegans), HvSKIP (barely) and OsSKIP (rice) had 
34, 13, 5, 1 and 35 interacted proteins, respectively [16,20- 
22,40,41]. Nevertheless, few SIPs could match each other 
among species. For example, both HvSKIP and GmGBPl 
interacted with GAMYB but OsGAMYB was not included 
in the 35 OsSIPs [21,22]. All the data indicated that SKIP 
might participate in distinct functions through the inter- 
action with diverse proteins. 

Conclusions 

In this study, GmGBPl ,as a homolog of SKIP in soybean, 
not only regulates plant flowering time but also alters 
plant resistance to abiotic stress in Arabidopsis. Such a 
flowering study about SKIP homolog in plant has no re- 
port so far. Although both OsSKIP and AtSKIP could 
improve plant tolerance to high salinity and drought 
[22,23], GmGBPl might improve plant tolerance to heat 
and drought, but reduce the resistance to high salinity. 
The flowering and stress related functions of GmGBPl 
might be used to develop soybean cultivars with the 
adaptability to broad grow regions and environment. 

Methods 

Plant material, growth conditions and chemical 
treatments 

Pure seeds of soybean cultivar 'DongNong 42' (photo- 
period-sensitive) were obtained from Soybean Research 
Institute of Northeast Agricultural University (Harbin, 
China). To analysis the tissue-specific expression of 
GmGBPl, the root, stem, leaf, trifoliate leaves, flower bud, 
pod and immature seed from soybean under both LDs 
and SDs were sampled. To analysis the expression of 
GmGBPl under abiotic stress conditions, the soybean 
plants, with the first fully expanding trifoliate leaves, were 
transplanted into 1/4 Murashige and Skoog (MS) [42] li- 
quid medium (no agar supplement) supplemented with 
100 [imol gibberellin, 100 \imo\ ABA, 200 mM NaCl or 
8% PEG6000. The plants that were transplanted into 1/4 
MS liquid medium without the supplements of agar, gib- 
berellin, ABA or NaCl were used as the control. The 
plants were sampled at 0, 0.5, 1, 2, 4 and 12 hours after 
transfer. 

Arabidopsis thaliana ecotype Columbia (Col-0) was 
used in this study. Seeds of atmyb33 T-DNA insertion 
line (Salk_058312) and atskip T-DNA insertion line 
(SAIL_681_H11) were obtained from The Arabidopsis 
Information Resource (TAIR). 

Identification of mutant 

PCR with the gene-specific primer pair and T-DNA- 
specific primer were used for genotyping to identify 
homozygous T-DNA inserted plants. The left genomic 



primer (LP), right genomic primer (RP) and the left T- 
DNA border primer (LB) for atskip were as follows: LP: 5 f 
CAAGCACAAGAGAGTCCCAAG 3' RP: 5' CGCCA 
CTTGCTCTCATAGTTC 3' LB: 5' GCCTTTTCA 
GAAATGGATAAATAGCCTTGCTTCC 3'. The LP, RP 
and LB primers for myb33 were as follows: LP: 5' ATCCA 
GAACTGTCAGACGCTG 3' RP: 5' AATTGCGTAT 
TTGGTTGGATG 3' LB: 5' ATTTTGCCGATTTCG 
GAAC 3 f . After confirmation of homozygous T-DNA in- 
sertion, gene knock-out was confirmed by RT-PCR with 
gene-specific primer. 

Gene expression analyses 

Total RNA was isolated from soybean leaves using TRIzol 
reagent (Invitrogen, USA). Real-time quantitative RT-PCR 
was performed on a Chromo4 Real-Time PCR System 
(Bio-Rad, USA) using SYBR Green PCR Master Mix 
Reagent (Takara, Japan). The measured Q values were 
converted to relative expression level using the AAQ 
method. Soybean Actin4 gene was used as the endogenous 
control. Primers used were qGmGBPl-F and qGmGBPl-R; 
qGmACTIN4-F and qGmACTIN4-R (Table 1). 

Semi-quantitative RT-PCR was used to analyze the ex- 
pression of flowering-related genes. Arabidopsis 18S 
rRNA gene was used as the endogenous control. The 
PCR products were analyzed through agarose gel elec- 
trophoresis and stained with ethidium bromide. All RT- 
PCR analyses were repeated three times, and one repre- 
sentative was shown. 

Yeast two-hybrid assay 

PCR-amplified cDNA fragments of GmGBPl, GmGBPla 
(amino acids 1-189), GmGBPlb (amino acids 190-356) 
and GmGBPlc (amino acids 357-612) were cloned 
into EcoRl and Pstl site of pBD- GAL4 as the bait con- 
structs; PCR-amplified cDNA fragments of GmGAMYBl, 
GmGAMYBla (amino acids 1-141) and GmGAMYBlb 
(amino acids 142-538) were cloned into BamHl and Pstl 
site of p AD -GAL4 as the prey constructs. The primers 
used were: yGmGBPl-F and yGmGBPl-R; yGmGBPl^-Y 
and yGmGBPla-R; yGmGBPlb-Y and y GmGBPl b-R; 
yGmGBPlc-Y and yGmGBPlc-R; yGmGAMYBl-Y and 
yGmGAMYBl-R; yGmGAMYBla-Y and yGmGAMYBla- 
R; yGmGAMYBlb-Y and yGmGAMYBlb-R (Table 1). All 
the bait constructs were transformed into the yeast strain 
YRG-2 respectively to detect the transcriptional activities. 

The yeast two-hybrid assay was performed using the 
GAL4 Two-Hybrid phagemid vector System (Stratagene, 
USA). The yeast strain YRG-2 was transformed with the 
bait plasmid, and the cells containing the bait constructs 
that have no self-activation activity were transformed 
with the prey plasmid according to the manufacturer s 
instruction manual. The co-transformants were selected 
on synthetic complete selection medium lacking Leu, 
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Table 1 Primers used in this study 



Primer names Primer sequences (5 -3 ) 



cGmGBP 7 -F 5'-TCCCAAAATCACCATCTA-3' 

cGmGBP 7 -R 5'-TOTAAC ATCC ATAAGCAGT-3' 

GmGBP 7-i-F 5'-GGGGACAAGmGTACAAAAAAGCAGGCTGCTATGGG 
GAAGAGGAGTAA-3' 

GmGBP 7-i-R 5'-GGGGACCACmGTACAAGAAAGCTGGGTCA™CCA 
TCCCGCATAG-3' 

yGmGBPl-V 5'-GAATOATGGCCACTCTGAAAGAGOTOTC-3' 

yGmGBPl-R 5'-GCGGCCGCCTAATGCCCTCmCAAATCCAATG-3' 

yGmGBPla-V 5'-GAATOATGGCCACTCTGAAAGAGOTOTC-3' 

yGmGBPla-R 5'-GCGGCCGCCTAATGCCCTCmCAAATCCAATG-3' 

yGmGBPl b-F 5'-GAATCGCGAAGTATATAAAGTACAAACCC-3' 

yGmGBPl b-R 5'-GCGGCCGCAATOTCTCCCCTCCAATOmCA-3' 

yGmG£P7c-F 5'-GAATCGGGGTOTACCAGCCGCGCCACCAG-3' 

yGmG8P7c-R 5'-GCGGCCGCCTAATGCCCTCmCAAATCCAATG-3' 

yGmGAMYB 1-F 5'-GGATCCATGAGACGAATGAAGAAAGATATO-3' 

yGmGAMYB 7-R 5'-GAATC™GAGGGGCTGGAATGGA^CA-3' 

yGmGAMYB 7 a-F 5'-GGATCCATGAGACGAATGAAGAAAGATATO-3' 

yGmGAMYB 7 a-R 5'-GAATCC^GATCCGGGTGTCCAGTAGT-3' 

yGmGAMYB 7 b-F 5'-GGATCCAACGGGCTGGGTOCCACmATCC-3' 

yGmGAMYB 7 b-R 5'-GAATO™GAGGGGCTGGAATGGA^CA-3' 

qGmGBPl-V 5'-mGTGAAGGAGAGTAGGGAGGAG-3' 

qGmGBPl-R 5'-™GTAGAGGCCATACCAAGAGCA-3' 

qGm/\C77/V4-F 5'-GTGTCAGCCATACTGTCCCCAm-3' 

qGm/\C77/V4-R 5'-GmCAAGCTOTGCTCGTAATCA-3' 

/\fS/C/P-F 5'-GGCTGmCAATGCGTOC-3' 

/US/C/P-R 5'-TOGGTCAACACmCACTCCTA-3' 

AtLFY-F 5 '-TGTG AACATCGC1TGTCGTC-3' 

AtLFY-R 5'-TAATACCGCCAACTAAAGCC-3' 

/\tC0-F 5'-AAGGTGATAAGGATGCCAAGGAG-3' 

/\fC0-R 5'-GGAGCCATAmGATATOAACTGA-3' 

/\f/T-F 5'-TGGTGGAGAAGACCTCAGGAAC-3' 

/\t/T-R 5'-TGCCAAGCTGTCGAAACAATAT-3' 

/\tS0C7-F 5'-CGTAAACTOTGGGAGAA-3' 

/\tS0C7-R 5'-TCAGAAOTGGGCTACTC-3' 

/\tF/_C-F 5'-TOTCGTCGTCTCCGCCTCC-3' 

AtFLC-R 5'-TOCTCCAGTOAACAAGAGCATC-3' 

/\tSl/P-F 5'-CCACATCCACCGACCATC-3' 

AtSVP-R 5'-AGGACOTGCCAAGCAGT-3' 

/\tMV / 833-F 5'-AGCGGAAACCTCTOTGGAT-3' 

/\tMVB33-R 5'-TCAmCCTATGTCATCGCC-3' 

AtFKF-F 5'-CAGCCATOGATGATACm-3' 

AtFKF-R 5'-GAGACAACCAGmAGAGCC-3' 

AtZAJ7-F 5 '-TGTCCG ATATGTGG AGTG A-3' 

AtZAT7-R 5'-TGTAACCAACGAGCCTGA-3' 

AtZAT12-F 5'-GCGAGTCACAAGAAGCCTAA-3' 



Table 1 Primers used in this study (Continued) 

AtZAT12-R 5'-ATCGGAAACTCCACTCCACAT-3' 

AtHSFA 7 6-F 5'-GmGGCTCATCATCCTC-3' 

AtHSFA 7 6-R 5'-TCCCTCCAGTAACTCATOA-3' 

AtHSFA2-F 5'-AAGGCGTOAACAATCCG-3' 

AtHSFA2-R 5'-CAGCGAACAACAmCCATA-3' 

AtRD26-F 5'-AATGGGTCGTCATCGTCT-3' 

/\f/?D26-R 5'-GCATCGTAACCACCGTAA-3' 

/AfMVC2-F 5'-GAACGACCCGTCTATGTGG-3' 

AtMYC2-R 5'-mCGG™TOTGOTGA-3' 

AtMYB96-F 5'-CmGOTGACCGTO^-3' 

AtMYB96-R 5'-CCTGATGGGTATGGA™-3' 

AtMYB2-F 5'-AAAAGCAAGCCAAACACC-3' 

AtMYB2-R 5'-CACTAATCTCGGCATCCA-3' 

At18SrRNA-F 5'-CGTCCCTGCCCmGTACAC-3' 

Atl8SrRNA-R 5'-CGAACAOTCACCGGATCA^-3' 



Trp and His. Large yeast clones appearing within 7 days 
were picked out to test the LacZ reporter gene. Positive 
clones were identified by the restriction digest analyses. 



Gene constructs and Arabidopsis transformation 

The GmGBP 1 -overexpression construct was constructed 
by directionally inserting the full cDNA into the cloned 
vector pMD18T and then into the binary vector pCAM- 
BIA3300. To make a dsRNAi construct of GmGBP 1, a 
505 bp fragment of GmGBP 7, 73% of nucleotide hom- 
ology to AtSKIP, was generated by PCR with primers 
GmGBPl-i-F and GmGBP 1 A-R (Table 1) and cloned 
into pJawohl8 vector through attB x attP (BP) recombin- 
ation cloning. The attBl and attB2 were the two sequences 
for the BP recombination reaction (Invitrogen, USA). The 
constructs were transformed into the Agrobacterium tume- 
faciens strain LBA4404 and used to transform Arabidopsis 
ecotype Columbia using the floral dip method (Clough and 
Bent 1998). 



Flowering analysis 

Seedlings were germinated on MS medium, and then trans- 
ferred to 1:1 of vermiculite and turfy-soil and grown under 
LDs or SDs conditions. The flowering time of Arabidopsis 
was measured by scoring the number of rosette leaves at 
the flowering time and the number of days from germin- 
ation to bolting [31]. At least 20 plants were analyzed each 
time, and the experiments were repeated for three times. 
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